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Abstract 

The dimensionally regularized massless on-shell planar triple box Feyn- 
man diagram with powers of propagators equal to one is analytically 
evaluated for general values of the Mandelstam variables s and t in a 
Laurent expansion in the parameter e = (4 — d)/2 of dimensional reg- 
ularization up to a finite part. An explicit result is expressed in terms 
of harmonic polylogarithms, with parameters and 1, up to the sixth 
order. The evaluation is based on the method of Feynman parameters 
and multiple Mellin-Barnes representation. The same technique can 
be quite similarly applied to planar triple boxes with any numerators 
and integer powers of the propagators. 
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In the last four years, the problem of analytical evaluation was completely solved 
for most important classes of two-loop Feynman diagrams with four external lines 
within dimensional regularization [lj. In the pure massless case with all end-points 
on-shell, i.e. p 2 = 0, i = 1, 2, 3, 4, this was done in El HI El EH U\ ■ The corresponding 
analytical algorithms were successfully applied to the evaluation of various two-loop 
virtual corrections jSJ. In the case of massless two- loop four-point diagrams with 
one leg off-shell the problem of evaluation was solved in [J0|, with subsequent 
applications [TT] to the process e + e~ — > 3jets. A first result for the massive on-shell 
case was presented in ^2]- (See ^Sl Ej for brief reviews of results on the analytical 
evaluation of various double-box Feynman integrals and the corresponding methods 
of evaluation.) 

In [03 E] , first analytical results on three-loop on-shell massless four-point dia- 
grams within dimensional regularization were obtained. The leading power asymp- 
totic behaviour of the dimensionally regularized massless on-shell planar triple box 
diagram shown in Fig. 1 in the Regge limit t/s — > was analytically evaluated in 
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Figure V. Planar triple box diagram. 

[To] with the help of the strategy of expansion by regions ^B]. Then, in |TJ], explicit 
analytical results for the unexpanded master planar triple box were presented for 1 / e J 
terms of Laurent expansion in e with j=6,5,4,3 and 2. 

The purpose of this paper is to complete this task i.e. analytically evaluate the 
missing 1/e part and the finite part. An explicit result will be expressed in terms of 
harmonic polylogarithms (HPL) with parameters and 1, up to the sixth order. 
The evaluation is based on the technique of alpha parameters and Mellin-Barnes 
(MB) representation which was successfully used in IU 13 U2] and reduces, due to 
taking residues and shifting contours, to a decomposition of a given MB integral into 
pieces where a Laurent expansion of the integrand in e becomes possible. 

The general planar triple box Feynman integral without numerator takes the form 



T(a 1 ,...,ai ;s,t;e) 
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where s = (pi + p 2 ) 2 and t — (p 2 + P3) 2 are Mandelstam variables and k, I, r are the 
loop momenta. Usual prescriptions k 2 = k 2 + iO, s = s + iO, etc. are implied. 

To resolve the singularity structure of Feynman integrals in e it is very useful to 
apply the MB representation 



/+ioo 
dz— -T(u + z)T(-z) , (2) 
-too y\ 



(X + Y) v T(v) 2ffi 



that makes it possible to replace sums of terms raised to some power by their products 
in some powers, at the cost of introducing extra integrations. By a straightforward 
generalization of two- loop manipulations Ej) one can introduce, in a suitable 
way, MB integrations, first, after the integration over one of the loop momenta, r, 
then after the integration over /, and complete this procedure after integration over 
the loop momentum k. As a result, one arrives 1 •"> at the following sevenfold MB 
representation of (JTJ): 

(m d / 2 ) 3 (-l) a 

T(o ll ... lfl8 ;M;6) = n j=2 ,5,7,8,9,ior(a,)r(4 - a 589{10 ) - 2e)(-s)-™* 

1 p dw ,ty r ( a 2 + W n- W )T( Z2 + Zi) T( Z3 + zA 

(2ffiy J-ioo jJ- 2 \sj r(ai + z 3 + z 4 )r(a 3 + z 2 + z±) 

T(2 - ai2 - e + z 2 )T(2 - a 23 - e + ^)r(a 7 + w - z 4 )r(-2: 5 )r(-z 6 ) 
r(4 - a 123 - 2e + iu - z A )T(a 6 - z 5 )r(a 4 - z 6 ) 
xT(+ai23 - 2 + e + z A )T(w + z 2 + z 3 + z A - z 7 ) 
xT(2 - a 5 9(io) -e- z 5 - z 7 )T(2 - a 589 - e - z 6 - z 7 ) 
xr(a 467 -2 + e + w-z 4 -z 5 -z 6 - z 7 )T(a 9 + -2 7 )r(a 5 + z b + z 6 + z 7 ) 
xT(4 - a 467 - 2e + z b + z 6 + z 7 )T(a bm{l0) - 2 + e + z 5 + z 6 + z 7 ) 
xT(2 - a e7 - e - w - z 2 + z 5 + z 7 )T(2 - a i7 - e - w - z 3 + z 6 + z 7 ) , (3) 

where a = Z^i^i) a 58 9(i ) = a 5 + a s + a 9 + a w , a 123 = ax + a 2 + a 3 , etc., and integration 
contours are chosen in the standard way. 

In the case of the master triple box, we set dj = 1 for i — 1, 2, . . . , 10 to obtain 



T^{s,t;e) = T(l,...,l;s,t;e) 



(iff 



d/2y 
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xT(l + e + w - z 4 - z 5 - z e - z 7 )T(—e - w - z 2 + z 5 + z 7 ) 

xT(-e -w- z 3 + z 6 + z 7 )T(1 + z 5 + z 6 + z 7 )T(2 + e + z 5 + z 6 + z 7 ) . (4) 

Observe that, because of the presence of the factor r(— 2e) in the denominator, we 
are forced to take some residue in order to arrive at a non-zero result at e = 0, so 
that the integral is effectively sixfold. 

Then the standard procedure of taking residues and shifting contours is applied, 
with the goal to obtain a sum of integrals where one may expand integrands in 
Laurent series in e. One- and two-loop examples of such procedures can be found, 
e.g., in ^31- The poles in e are not visible at once, at a first integration over one 
of the MB variables. However, the rule for finding a mechanism of the generation 
of poles is based on the simple observation that a product of two gamma functions 
Via + z)T(b — z), where z is a MB integration variable and a and b depend on the 
rest of the variables, generates a pole of the type T(a + b). This means that any 
contour in the next integrations should be chosen according to this dependence. So, 
the first step is an analysis of various pairs of gamma functions and various orders 
of integration in (@J. The analysis of the integrand shows that the following four 
gamma functions play a crucial role for the generation of poles in e: T(— e + ^2,3) 
and T(— 1 — e — zq^ — z 7 ). The first decomposition of the integral (jlj) arises when 
one either takes a residue at the first pole of one of these gamma functions or shifts 
the corresponding contour, i.e. changes the nature of this pole. As a result (@J) is 
decomposed as 2T 000 i + 2T 0010 + 2T 00 n + T 010 i + 2T OU0 + 2T 0111 + T 10 i + 2T lon + T im 
where a symmetry of the integrand is taken into account. Here the value 1 of an 
index means that a residue is taken and means a shifting of a contour. The first 
two indices correspond to the gamma functions T(— e + z 2 ) and T(— 1 — e — z 5 — z 7 ) 
and the second two indices to T(— e + Z3) and T(— 1 — e — zq — z 7 ), respectively. The 
term T ooo is absent because it is zero at e = due to T(— 2e) in the denominator. 

Each of these terms is further appropriately decomposed and, eventually, one is left 
with integrals where integrands can be expanded in e. These resulting terms involve 
up to five integrations. Taking some of these integrations with the help of the first 
and the second Barnes lemmas, one reduces all the integrals to no more than twofold 
MB integrals of gamma functions and their derivatives. In some of them, one more 
integration can be also performed in gamma functions. Then the last integration, 
over w is performed by taking residues and summing up resulting series, in terms 
of HPL. Keeping in mind the Regge limit, t/s 0, let us, for definiteness, decide 
to close the contour of the final integration, over w, to the right and obtain power 
series in t/s. The coefficients of these series are (up to (— l) n ) linear combinations of 
1/n 6 , S\(n)/n 5 , . . . , S\{n)S 3 {n) /n 2 , . . ., where Sfc(n) = YIj=i3~ k - Summing up these 
series gives results in terms of HPL of the variable —t/s which can be analytically 
continued to any domain from the region \t/s\ < 1. 

In the twofold MB integrals where one more integration (over a variable differ- 
ent from w ) can be hardly performed in gamma functions, one performs it with w 
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in a vicinity of an integer point w = n = 0, 1, 2, . . ., in expansion in z = w — n, 
with a sufficient accuracy. Then one obtains powers series where, in addition to 
l/n Q ,Si(n)/n b ,..., quantities like S ik (n) = J2]=ij~ l S k (j), S ik i(n) = T,] =1 j~ l S k i(j) 
appear. These series are also summed up in terms of HPL. 
Eventually we arrive at the following result: 

where 7e is the Euler constant, x = —t/s, L = ln(s/t), and 

16 5 r 3 2 

ce = y, c 5 = --L, c 4 = --7T , 

3 11 1 31 

c 3 = 3(tf ,o,i(x) + LH 0A (x)) + -(L 2 + n^H.ix) - -ix 2 L - — ( 3 , 

c 2 = -3 (17iJ ,o,o,i(x) + #o,o,i,i(x) + # ,i,o,i(x) + #1,0,0,1 ( x )) 
-L (37# ,o,i(x) + 3^0,1,1^) + 3#i,o,i(x)) " \(L 2 + 7i 2 )H hl (x) 

- (f L 2 + 8tt 2 ) H 0A (x) - (Jl 3 + tt 2 L - 3C 3 ) #i(x) + f Cs^ - ^ 4 , 

ci = 3 (81if ,o,o,o,i(x) + 41 #o,o,o,i,i(x) + 37 - H o,o ) i,o,i( a: ) + #0,0,1,1,1 ( x ) 
+33# ,1,0,0, i(x) + #0,1,0,1, i(x) + #o,i,i,o,i(x) + 29^1^,0,0,1 (x) 
+#1,0,0,1,1 (x) + #i,o,i,o,i(x) + #i,i,o,o,i(x)) + L (177.fZo,o,o,i(x) + 85^0,0,1,1 (x) 
+73^0,1,0,1 (x) + 3if ,i,i,i(x) + 61i7i, ,o,i(x) + 3iJi, ,i,i(x) + 3#i,i, ,i(x)) 
+ ^119 L2 + 139^ ^ ^ } + ^ + ^ 

+ (|L 2 + 14tt 2 ) H lflA (x) + l(L 2 + tt 2 ) ^ 1>lfl (x) 

+ (f L 3 + ^n 2 L - 96C 3 ) Ho,i(x) + (Jl 3 + n 2 L - 3Cs) #i,i(x) 

/9 r4 25 2 r2 roA r 13 A rr . , 503 4r 73 2a 301 a 

+ I -L + — 7T 2 L 2 - 58C 3 £ + — 7T 4 if^x) 7T 4 L + — 7T 2 C 3 Cs , 

V8 8 S 8 J W 1440 4 S3 15 S ' 

c = - (951if ,o,o,o,o,i(x) + 819# ,o,o,o,i, i(x) + 699if ,o,o,i,o,i(x) + 195if ,o,o,i,i,i(x) 

+547iJ ,o,i,o,o,i(x) + 231iJ ,o,i,o,i,i(x) + 159if ,o,i,i,o,i(x) + 3if ,o,i,i,i,i(x) 

+363iJ ,i,o,o,o,i(x) + 267iJ ,i,o,o,i,i(x) + 195if ,i,o,i,o,i(x) + 3if ,i,o,i,i,i(x) 

+123if ,i,i,o,o,i(x) + 3i? ,i, 1,0,1, i(x) + 3i? ,i,i,i,o,i(x) + 147-ffi, , ,o,o,i(x) 

+303iJi, ,o,o,i,i(x) + 231ifi )0 ,o,i,o,i(x) + 3 #1,0,0,1, i,i(x) + 159-ffi, ,i,o,o,i(x) 

+3#i, ,i,o,i,i(x) + 3^i,o,i, 1,0,1 (a?) + 87iJi,i, , ,o,i (x) + 3iJi,i, , ,i,i(x) 

+3ifi,i, ,i,o,i(x) + 3i?i,i,i, ,o,i(x)) 

-L (729# ,o,o,o,i(x) + 537# ,o,o,i,i(x) + 445# ,o,i,o,i (x) + 133# ,o,i,i,i(x) 
+321# ,i,o,o,i(x) + 169# ,i,o,i,i(x) + 97# ,i,i,o,i(x) + 3# ,i,i,i,i(x) 
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+165# 1A o, 0) i(V) + 205^1,0,0,1,1 (x) + 133#i A i A iO) + 3#i, ,i,i,i(z) 
+61^1,1,0,0,1 (x) + 3^1,1,0,1,1(0;) + 3^1,1,1,0,1(2;)) 

- (^L 2 + ^vr 2 ) H 0fiA1 {x) - (^L 2 + ^vr 2 ) H 0>0)hl (x) 

- (jj-L* + ^ 2 ) ff W (i) - (f^ 2 + 32vr 2 ) 

- - ^vr 2 ) Fi, , ,i(x) - (±^L 2 + sovr 2 ) F W|1 (x) 



- (^yL 2 + 147T 2 j ifl,l,o,l(x) - - [L 2 + 7T 2 j #l,l,l,l(x) 

- + ^ 2 £ - 455C 3 ) H ofiA (x) - (yL 3 + ^vr 2 L - 120C 3 ) H 0>1>1 (x) 
f £ 3 + - 156C 3 ) F 1A i(x) - (Jl 3 + tt 2 L - 3C 3 ) #i,i,i(*) 
f + ^vr 2 L 2 - 291(3^ + ^vr 4 ) ^(s) 
^ 4 + fvr 2 L 2 -58C3L + fvr 4 )i/i,i(x) 

27 r5 25 2r3 183 ^ r2 131 4 r 37 2 ^ rr7A \ rr , , 

— L H vr 2 L 3 (3^ H ttL 7T %3 + 57C 5 #i(x) 

40 8 2 S3 60 12 S3 S5 / u ; 

/223 2> , \ 167 a2 624607 fi 

+ tt 2 C 3 + 149C 5 L + Cl . 6 

V 12 S3 S V 9 S3 544320 V ; 

Here (3 = C(3),Cs — C(5) and ((z) is the Riemann zeta function. The functions 
H ai ,a 2 ,...,a n (x) = H(ax, 02, ... , o n ; sc), with Oj = 1, 0, —1, are HPL [T2j which are recur- 
sively defined by 

rx 

H(a u a 2 , . . . , a n ; x) = /(ai; t)H(a 2 , ...,a n ;t), 

Jo 

where 

/(±l;x) = 1/(1 =fx), f(0;x) = l/x, H(±l;x) = tH 1 T x), H(0; x) = \nx . 

In (pj). only HPL with parameters and 1 are involved. If a given HPL involves 
only parameters Oj = and 1 and the number of these parameters is less or equal 
to four, it can be expressed [T7j in terms of usual polylogarithms Li (a;) [THj and 
generalized polylogarithms 

s. t w= ( r ir t: /"-- 1(t) '" t(i -*v 



(o- 1)!6! Jo t 

(See JH] where the coefficients cy, with j > 2 are expressed in terms of (generalized) 
polylogarithms.) 
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The above result was confirmed [201 with the help of numerical integration in the 
space of alpha parameters [2"T l. Another natural check of the result is its agreement 
with the leading power Regge asymptotic behaviour ^Hj which was evaluated by an 
independent method based on the strategy of expansion by regions [TH] . 

The procedure described above can be applied, in a similar way, to the calculation 
of any massless planar on-shell triple box. At a first step, one has to take care of the 
following four gamma functions in (|3*J): 

T(2-a 1 2-e + Z2),T(2-a 2 3-e + z 3 ),T(2-a 59(w) -e-z 5 -z 7 ),T(2-a 589 -e-ZQ-z 7 ) . 

This procedure gives a decomposition similar to 2T 000 i + 2T 00 io + — Next steps will 
be also generalizations of the corresponding steps in the evaluation of (JU). Hopefully, 
such a procedure can be made automatic by means of computer algebra. 

The result presented above shows that analytical calculations of four-point on-shell 
massless Feynman diagrams at the three-loop level are quite possible so that one may 
think of evaluating three-loop virtual corrections to various scattering processes. 
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